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1.0  INTRODUCTION 

The  work  reported  herein  was  performed  by  the  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Materiel  Command  (AFMC),  under  Program  Element  65807F.  The  Air  Force 
Program  Managers  were  Capt.  D.  G.  Burgess,  Maj.  H.  Martin,  and  Capt.  P.  Zeman,  DOT.  The  work 
was  performed  by  Micro  Craft  Technology,  support  contractor  for  aerodynamic  testing  at  AEDC, 
AFMC,  Arnold  Air  Force  Base,  TN.  The  work  was  performed  in  the  Aerospace  Systems  Facility 
(ASF)  and  the  Technology  and  Development  Facility  (TDF)  under  AEDC  Project  Number  DDOl 
(Job  Number  0115).  The  work  was  conducted  during  the  period  between  1  October  1989  and  30 
September  1994. 

Electric  arc  heaters  have  been  used  at  the  AEDC  for  the  simulation  of  reentry  flight  heating  at 
pressures  up  to  150  atm  and  mass  average  enthalpy  of  4.65  -  9.3  x  10®  J/kg  (2000  -  4000  Btu/lbm)  in 
air.  Corresponding  heat  flux  levels  to  the  wall  of  the  arc  heater  nozzles  are  as  high  as  9  x  10’  W/m’ 
(8,(XX)  Btu/fF  sec),  and  future  test  conditions  will  require  nozzles  to  survive  2(X)  atm  (Ref.  1)  at  heat 
flux  levels  up  to  1.25  x  10*W/m^(ll,0(X)Btu/fF  sec).  Backside  subcooled,  forced  convection  cooling 
with  high-pressure  demineralized  water  is  currently  used  to  balance  the  heat  load  originating  from  the 
high-temperature  gas  flowing  through  the  nozzle  (Fig.  I).  In  order  to  accommodate  the  higher  heat 
loads,  there  is  current  interest  in  exploiting  the  backside  water  cooling  concept  to  the  limit  of  its 
capability.  In  addition  to  air-side  heating  prediction  and  favorable  material  characteristics,  the 
thermostructural  design  of  advanced  nozzle  concepts  requires  reasonable  estimates  of  cooling  heat- 
transfer  coefficients  and  limitations  (i.e.,  burnout  limit  or  critical  heat  flux  condition)  of  water  cooling, 
as  discussed  in  Ref.  2. 

Efficient  cooling  of  arc  heater  nozzles  is  accomplished  with  high  heat-transfer  coefficients 
resulting  from  the  transition  of  the  coolant  from  single-phase  convection  to  nucleate  boiling  at  high¬ 
er  heat  flux.  Analytical  modeling  of  the  heat-transfer  mechanisms  for  the  design  of  complex 
configurations  becomes  difficult  in  the  boiling  regime,  especially  at  the  critical  heat  flux  (CHF) 
condition.  Experimental  investigation  of  the  cooling  process  can  provide  the  means  to  study  the 
heat-transfer  mechanisms,  evaluate  parametric  trends,  and  develop  working  correlations  for  the 
cooling  configuration  of  interest.  A  flow  boiling  apparatus,  called  the  High  Temperature  Wall 
Laboratory  (HTWL),  has  been  developed  at  the  USAF/AEDC  to  perform  experimental 
investigation  of  the  cooling  processes  encountered  in  high-pressure,  electric  arc  heater  facilities.  A 
summary  of  the  development  and  operation  of  the  apparatus  and  a  discussion  of  initial  experimental 
work  using  the  apparatus  are  contained  in  this  report. 
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2.0  BACKGROUND 

Once  a  cooling  process  has  transitioned  into  boiling,  theoretical  approaches  for  determining 
heat-transfer  characteristics  have  limited  application  because  of  difficulties  in  obtaining  interaction 
or  interface  properties  of  the  two  phases.  Moreover,  no  theory  has  yet  been  created  for  forced 
convection,  or  flow,  boiling  burnout  (Ref.  3).  Many  researchers  have  resorted  to  developing 
correlations  based  on  microscopic  bubble  mechanisms,  dimensional  analysis,  fairing  of  test  data,  or 
a  combination  of  these.  Reasonable  confidence  in  predicting  heat-transfer  information  in  the  partial 
nucleate  boiling  and  fully  developed  nucleate  boiling  regimes  has  been  shown  by  Bergles  and 
Rohsenow,  Ref.  4,  and  Rohsenow,  Ref.  5,  respectively  (see  also  Guglielmini,  et  al,  Ref.  6); 
however,  considerable  disagreement  between  the  large  number  of  CHF  correlations  exists.  Gambill 
(Ref.  7)  compared  several  CHF  correlations  for  water  flow  in  a  tube  and  found  that  prediction  of 
CHF  at  higher  coolant  velocity,  using  the  correlations,  varied  by  nearly  a  factor  of  four  and  by  a 
factor  of  two  at  lower  velocity.  Factors  of  two  or  greater  between  correlation  predictions  at  higher 
coolant  mass  velocity  have  also  been  noted  by  Zeigamik,  et  al  (Ref.  8)  and  Boyd  (Ref.  9),  and  the 
disagreement  tends  to  worsen  as  mass  velocity  increases.  Empirical  CHF  correlations  typically 
have  been  derived  for  a  specific  range  of  experimental  conditions,  and  extrapolation  of  a  correlation 
outside  the  specific  range  can  lead  to  very  large  errors.  In  addition,  a  given  analytical  CHF  correla¬ 
tion  typically  has  enough  “adjustable”  constants  in  the  predictive  equations  to  permit  an  acceptable 
correlation  of  a  specific  data  set;  however,  general  application  of  the  correlation  to  various  config¬ 
urations  or  test  conditions  may  result  in  significant  errors  (Ref.  10).  As  recommended  in  Ref.  11, 
among  the  more  appropriate  CHF  correlations  for  subcooled,  forced  convection  water  flow  are 
those  proposed  by  Bemath  (Ref.  12),  Van  Huff  and  Rousar  (Ref.  13),  Rousar  (Ref.  14),  Yagov  and 
Puzin  (Ref.  15),  Levy  (Ref.  16),  and  Labuntsov  (Ref.  17). 

Various  flow  boiling  experiments  performed  in  the  past  have  provided  insight  to  the  boiling 
heat-transfer  mechanisms  and  parametric  trends  which  have  aided  in  the  development  of  working 
correlations.  Resistance  heating,  or  Joule  heating,  has  been  used  by  researchers  (boiling  heat 
transfer)  for  a  number  of  years  to  produce  the  necessary  heat  flux  from  a  surface  to  be  cooled  (Ref. 
11).  The  heat  flux  level  is  more  easily  determined  from  a  simple  energy  balance  in  a  Joule  heating 
configuration,  where  other  heating  techniques  such  as  secondary  fluid  heating  (e.g.,  steam)  or 
cartridge  heaters  introduce  more  complex  energy  conversion.  Though  numerous  external  flow  boil¬ 
ing  experiments  have  been  performed  (Ref.  3),  boiling  heat  transfer  is  sensitive  to  the  coolant  flow 
geometry  and  heated  wall  configuration;  therefore,  internal  flow  boiling  data  on  configurations 
closely  representing  arc  heater  nozzle  cooling  passages  are  more  appropriate  for  the  analysis  of 
cooling  processes  in  an  arc  heater  configuration.  Internal  flow  boiling  has  been  studied  primarily 
with  tubes  and  annuli,  where  significant  differences  in  heat-transfer  data  have  been  noted  between 
tubes  and  annuli  of  the  same  length  and  equivalent  diameter  (e.g..  Refs.  18  and  19).  In  addition, 
more  than  20  parameters  have  been  identified  that  affect  the  CHF  in  flow  boiling  (Ref.  20).  Further 
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flow  boiling  experimentation  is  required  because  of  the  sensitivity  of  boiling  heat-transfer  data  to 
geometry  and  test  conditions,  and  the  large  uncertainty  resulting  from  the  extrapolation  of  current 
boiling  heat-transfer  correlations  to  the  arc  heater  cooling  requirements. 

3.0  EXPERIMENTAL  APPARATUS 

In  order  to  verify  the  predictive  approaches  and,  if  necessary,  develop  a  more  appropriate 
CHF  correlation  for  the  arc  heater  nozzle  cooling  conditions,  an  experimental  boiling  apparatus, 
called  the  High  Temperature  Wall  Laboratory  (HTWL),  was  developed  at  the  AEDC.  The  desired 
heat  flux  at  the  surface/coolant  interface  on  the  inner  wall  of  an  annulus  is  achieved  by  passing  high 
current  through  a  thin-walled  tubular  metallic  test  section,  thereby  producing  Joule  heating.  The 
annular  configuration  was  selected  for  the  experimental  studies  since  it  represents  the  coolant  flow 
geometry  of  the  current  arc  heater  nozzles.  Heat  flux  levels  as  high  as  2  x  10*  W/m^  (18,000  Btu/ffl 
sec)  are  desired  in  the  apparatus  for  the  study  of  future  nozzle  concepts.  This  maximum  heat  flux 
level,  coolant  flow  requirements  anticipated  in  advanced  nozzle  concepts,  instrumentation  require¬ 
ments  for  the  definition  of  operating  conditions,  and  reasonable  machining  size  limitations  dictate 
the  overall  size  and  configuration  of  the  test  section  and  support  systems.  A  schematic  of  the 
apparatus  and  support  systems  is  shown  in  Fig.  2. 

3.1  POWER  REQUIREMENTS  AND  ELECTRICAL  SYSTEM 

The  peak  heat  flux  occurs  near  the  throat  of  an  arc  heater  nozzle  where,  for  the  baseline 
nozzle  configuration  under  consideration,  the  inner  diameter  of  the  annular  cooling  channel  is 
approximately  1  in.  (Ref.  1).  By  selecting  a  low-cost,  machinable,  low  electrical  conductance 
metal  such  as  304  stainless  steel  (material  properties  given  in  Appendix  A),  and  following  the 
computation  approach  outlined  in  Appendix  B,  the  power  supply  requirements  necessary  to 
achieve  up  to  2  x  10*  W/m^  (18,000  Btu/ff  sec)  in  a  25.4-mm  (1-in.),  19-mm  (0.75-in.),  and  a 
12.7-mm  (0.5-in.)  diam  tube  are  approximated,  as  shown  in  Fig.  3*.  As  pointed  out  in  Ref.  1 1,  the 
heated  wall  material  has  been  shown  to  have  a  possible  effect  on  boiling  heat  transfer;  therefore, 
consideration  of  the  actual  copper-zirconium  (Amzirc®)  nozzle  material  should  be  included.  Be¬ 
cause  copper  and  its  alloys  have  very  high  electrical  conductivity,  the  materials  require  very  high 
power  in  resistance  heating  to  reach  the  heat  flux  levels  identified  above.  Figure  4  illustrates  the 
high  power  requirements  for  a  25.4-mm  (1-in.),  19-mm  (0.75-in.),  and  a  12.7-mm  (0.5-in.)  diam 
Amzirc  tube. 


*This  analysis  assumes  a  mean  material  temperature  that  is  computed  from  the  HTWL  data  reduction  program 
discussed  in  Section  3.5,  Data  Reduction.  The  computation  of  the  mean  temperature  requires  an  energy  balance  of  the 
heated  wall  and  cooling  fluid,  which  is  accounted  for  in  the  heat-transfer  data  reduction.  This  mean  temperature  is  used 
here  to  determine  reasonable  values  of  the  material  resistivity.  Note  that  the  same  mean  temperature  is  used  for  the 
different  materials,  i.e.,  no  material  effects  on  the  heat-transfer  process  are  assumed. 
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Experimental  results  have  shown  that  an  a-c  type  of  power  source  can  adversely  affect  the 
boiling  process  (Ref.  11).  Therefore,  two  16,000-amp  d-c  rectifiers  manufactured  by  Rapid  Power 
Technologies  were  selected  for  the  HTWL  test  section  heating.  The  rectifiers,  when  operated  in  a 
parallel  mode,  would  provide  up  to  32,000  amps  at  100  volts  (3.2  MW)  for  heating  of  the  HTWL 
test  section.  Specifically,  the  high  heat  flux  condition  discussed  previously  could  be  achieved  with 
a  25.4-mm  (1.0-in.)  diam  304  stainless-steel  tube  with  a  0.76-mm  (0.030-in.)  thick  wall  (Fig.  3a). 
Material  effects  would  be  evaluated  with  12.7-mm  (0.5'in.)  diam,  0.25-mm  (0.010-in.)  thick  tubes^ 
of  various  materials  including  Amzirc  (Fig.  4c)  and  stainless  steel  (Fig.  3c).  One  disadvantage  of 
the  rectifiers  is  that  peak  performance  is  achieved  when  the  rectifiers  are  operated  at  the  peak  volt¬ 
age  output  of  100  volts.  The  test  section  length  to  be  heated,  material  type,  and  internal  temperature 
drive  the  voltage  drop  in  the  system;  therefore,  an  optimum  test  section  length  of  150  mm  (6.0  in.) 
was  selected,  and  an  adjustable  water-cooled  ballast  resistor  bank  was  added  to  the  electrical  system 
to  allow  operation  of  the  rectifiers  at  peak  voltage  output  for  improved  power  accuracy.  When  test 
sections  with  low  voltage  drop  are  used  (e.g.,  Amzirc  and  copper),  an  additional  voltage  drop  is  set 
at  the  ballast  resistor  bank  to  operate  the  rectifiers  near  100  volts.  The  actual  uncertainty  of  the  pow¬ 
er  levels  provided  by  the  rectifiers  is  discussed  in  Section  5.0,  Performance  and  Data  Uncertainty. 
Current  loss  into  the  cooling  water  of  the  resistor  bank  was  determined  to  be  negligible  (less  than  1 
mA  with  raw  water  cooling).  The  rectifiers  are  routinely  operated  up  to  100  percent  (1 .6  MW)  pow¬ 
er  in  the  current  control  mode.  Rectifier  output  from  1.6  MW  to  3.2  MW  requires  operation  in  the 
voltage  control  mode  with  limited  duration. 

The  mezzanine  of  the  High  Temperature  Laboratory  (HTL)  was  selected  as  the  site  for 
HTWL  primarily  because  of  the  proximity  of  the  6900-V  electrical  feed  required  by  the  rectifiers 
and  the  high-pressure  demineralized  water  system  in  the  building.  Figure  5  shows  the  layout  of  the 
HTWL  equipment  positioned  on  the  HTL  mezzanine.  The  6900-V  electrical  feed  provides  power 
to  the  rectifiers  through  a  pad-mounted  switch.  Power  for  the  rectifier  controls  is  supplied  by  the 
building  480- V  system.  The  high  current,  low  voltage  power  provided  by  the  rectifiers  is  fed  to  the 
ballast  resistor  bank  and  test  section  assembly  through  eight  38-mm  (1.5-in.)  diam,  777.7  MCM 
insulated  cables.  The  ballast  resistor  bank  is  made  up  of  eight  tubes  (four  for  each  rectifier)  made 
of  standard  19-mm  (0.75-in.)  diam  Sch-40  304  stainless-steel  pipe.  The  insulated  power  cables  are 
attached  to  the  tube  elements  with  copper  lugs.  Copper  jumper  plates  between  pairs  of  the  tubes 
allow  the  length  to  be  varied  from  zero  to  2.75  m  (9  ft).  Such  an  arrangement  allows  voltage  drops 
in  the  elements  from  zero  to  80  volts  to  be  set.  Each  pipe  is  electrically  isolated  from  ground  by 
Plexiglas®  supports  and  high-pressure  flexible  hoses  at  the  inlet  and  exit.  The  tube  elements  are 
cooled  by  low-pressure  (7  bar  or  100  psi)  raw  water  through  upstream  and  downstream  manifolds. 
Details  of  the  ballast  resistor  bank  are  shown  in  Fig.  6. 


^Extremely  thin  tubes  (less  than  0.25-mm  thick  wall)  are  difficult  to  instrument,  easily  damaged,  and  have  higher  volt¬ 
age  requirements  than  are  capable  with  the  selected  rectifiers  and  are,  therefore,  not  considered. 
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A  ground  fault  protection  system  prevents  the  rectifiers  from  being  started  should  a  ground 
fault  be  present  anywhere  in  the  high  energy  electrical  system.  In  addition,  the  ground  fault  system 
will  phase  back  and  shut  down  the  rectifiers  during  operation  if  a  ground  fault  occurs  at  any  power 
setting.  For  example,  during  initial  testing  of  a  large  diameter  test  section,  which  had  a  small  coolant 
flow  annulus,  the  test  section  failed  at  burnout  and  bent  against  the  stainless-steel  water  jacket. 
Because  the  test  section  shell  and  the  water  jacket  were  grounded  through  the  test  stand  (isolated 
later  with  Plexiglas  spacers),  a  ground  path  for  the  high  electrical  energy  occurred.  The  ground  fault 
protection  system  detected  the  fault  and  shut  down  the  rectifiers  before  significant  damage  to  the 
hardware  could  occur.  Because  small  current  paths  to  ground  typically  exist  in  the  HTWL  circuit 
(e.g.,  conductance  through  the  demineralized  water  adjacent  to  the  energized  test  section)  and  the 
fact  that  the  level  changes  with  rectifier  power  setting,  an  adjustable  trigger  point  for  the  ground 
fault  system  is  provided. 

3.2  COOLANT  REQUIREMENTS  AND  DEMINERALIZED  WATER  SYSTEM 

The  existing  closed-loop,  high-pressure  demineralized  water  system  used  by  the  AEDC  arc 
facilities  is  capable  of  high  volumetric  flows  (2500  liters  per  min  or  650  gpm)  at  a  pressure  up  to 
100  bar  (1500  psi).  Coolant  velocity  at  the  throat  of  the  baseline  nozzle  configuration  is 
approximately  30.5  m/s  (100  fps)  with  a  total  mass  flow  rate  of  approximately  5  kg/sec  (11  Ibm/ 
sec)  or  300  liters  per  min  (80  gpm).  Future  pumping  requirements  may  necessitate  higher  pressure 
coolant  flow;  however,  the  high  cost  for  such  a  system  in  a  continuous  flow  facility  will  probably 
limit  operation  in  the  foreseeable  future  to  coolant  pressure  below  100  bar.  Therefore,  the  closed- 
loop,  high-pressure  demineralized  water  system  is  adequate  for  constant  coolant  inlet  temperature 
studies  in  the  HTWL. 

The  closed-loop  circuit  (Fig.  2)  incorporates  an  existing  Bingham  high-pressure  pump  (Fig. 
7)  and  heat  exchanger  located  in  the  HTL,  and  connection  to  the  high-pressure  loop  for  HTWL 
supply/retum  is  below  the  mezzanine  near  the  H2  Arc  Heater  Facility.  Because  the  existing  closed- 
loop  water  system  pumps  approximately  2500  liters  per  min  (650  gpm),  a  bypass  loop  in  the  HTWL 
system  allows  mass  flow  variation  over  a  broad  range  with  the  velocity  being  set  at  the  test  section 
by  the  annular  flow  area.  Strainers  with  0.254-mm  (0.01-in.)  mesh  are  located  upstream  and  down¬ 
stream  of  the  test  section  to  prevent  contaminants  from  entering  the  test  section  or  returning  to  the 
pump.  Isolation  valves  in  the  HTWL/H2  flow  loops  permit  quick  changeover  between  the  two 
installations.  The  entire  circuit  is  constructed  of  304  stainless  steel  except  for  short  sections  of  high- 
pressure  flexible  hose  used  to  electrically  isolate  the  HTWL  test  section.  A  51 -mm  (2-in.)  diam 
throttling  valve  (Annin  Co.  globe  valve)  and  Grove  pressure  regulator  upstream  of  the  test  section 
and  a  51 -mm  (2-in.)  diam  flow  control  valve  (Kent  Introl  Ltd.  globe  valve  with  pneumatic 
positioner)  downstream  are  used  to  adjust  flow  conditions  at  the  test  section  and  maintain  flow 
stability  in  the  circuit  (see  discussion  in  Ref.  1 1  concerning  flow  loop  instabilities).  A  51-mm  (2- 
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in.)  diam  Hoffer  Flow  Controls  turbine  flowmeter  is  installed  upstream  of  the  coolant  pressure 
regulator  and  inlet  manifold.  High-pressure  pump  capabilities  permit  run  times  up  to  30  min.  The 
bulk  coolant  temperature  is  fixed  at  approximately  room  temperature,  and  the  water  conductivity  is 
limited  to  approximately  100  |iS/cm.  The  dissolved  oxygen  content  is  typically  4  ppm  or  31.5 
percent  of  saturation. 

3.3  TEST  SECTION  ASSEMBLY 

The  HTWL  test  section  assembly  (Fig.  8)  currently  in  use  simulates  the  annular  flow  arrange¬ 
ment  that  is  incorporated  in  the  cooling  passage  of  an  arc  heater  nozzle.  The  actual  arc  heater  nozzle 
and  associated  cooling  passages  are  converging/diverging  by  design  (Fig.  1),  and,  consequently, 
parameters  such  as  the  coolant  mass  velocity,  static  pressure,  subcooling,  acceleration,  and  heat  flux 
distribution  vary  along  the  length  of  the  nozzle.  Because  more  than  20  parameters  have  been  shown 
to  influence  flow  boiling  heat  transfer  and  CHF  (Refs.  1 1  and  20),  the  test  section  assembly  in 
HTWL  was  designed  such  that  individual  parametric  trends  could  be  evaluated.  Figure  8b  shows  a 
cutaway  of  the  test  section  assembly,  revealing  the  heated  horizontal  tubular  test  section  and 
annular  flow  channel. 

Electrical  energy  enters  the  assembly  by  way  of  eight  large  electrical  cables  clamped  to  four 
copper  tabs  on  one  end  of  the  test  section  assembly.  The  energy  passes  through  a  copper  flange  and 
test  section  end  piece  to  the  test  section  which  is  electron  beam  welded  to  the  copper  end  pieces. 
The  electrical  energy  passes  out  of  the  assembly  through  a  copper  flange  and  cable  attachments, 
similar  to  the  way  it  enters.  All  electrical  conducting  components  of  the  test  section  assembly, 
except  the  test  section  tube,  are  fabricated  from  electrolytic  tough-pitch  copper.  The  nonconducting 
portions  of  the  assembly  are  fabricated  from  304  stainless  steel  and  are  insulated  from  the  conduct¬ 
ing  components  by  Micarta®  and  C-1 1  glass/epoxy  insulating  materials.  In  addition,  the  entire  test 
section  assembly  is  insulated  from  the  test  stand/support  structure  by  Plexiglas  spacers.  Cool-Amp 
Silver  Plating  Powder  and  Bumdy  Penetrox  A  Anti-Oxidation  Coating  are  used  where  slip  fit/ 
electrical  contact  assembly  is  required.  As  shown  in  Fig.  8b,  the  downstream  test  section  end  piece 
includes  a  piston  assembly  to  allow  for  test  section  thermal  expansion  which,  in  some  cases,  can  be 
as  large  as  2.5  mm  (0.1  in.). 
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The  water  coolant  flow  enters  the  test  section  assembly  from  high-pressure  flexible  hoses  at  four 
fittings  spaced  between  the  electrical  connections  (see  Fig.  8b).  The  flow  channels  transition  into  a 
horizontal  annular  flow  configuration  that  includes  an  entry  length  of  at  least  ten  hydraulic  diameters 
for  fully  developing  the  flow  prior  to  reaching  the  heated  test  section.  Interchangeable  water  jacket 
sleeves  (an  example  is  shown  in  Fig.  9a)  allow  the  annulus  gap  to  be  varied.  Hydraulic  diameters^  of 
3.8  mm  to  7.1  mm  (0.15  in.  to  0.28  in.)  were  selected  based  on  the  test  section  sizes  discussed  in  the 
following  paragraph  and  a  velocity  range  requirement  of  1 5  to  61  m/s  (50  to  200  fps).  The  water  jacket 
sleeves  are  fabricated  of  304  stainless  steel  or  Delrin®  AF  resin/Teflon®  fiber  composite.  The  compos¬ 
ite  water  jacket  is  used  primarily  for  the  small  hydraulic  diameter  where  arcing  across  the  water  gap 
may  be  possible.  The  coolant  exits  the  test  section  assembly  similar  to  the  way  it  enters.  A  pressure 
relief  valve  connected  to  the  annular  flow  channel  protects  the  assembly  from  over-pressurization 
when  large  vapor  voids  are  generated  at  test  section  burnout. 

Three  sizes  of  test  section  tubes  are  currently  in  use:  a  26-mm  (1 .024-in.)  diam  tube  with  a  0.76- 
mm  (0.03-in.)  thick  wall,  a  19-mm  (0.75-in.)  diam  tube  with  a  0.51-mm  (0.02-in.)  thick  wall,  and  a 
12-mm  (0.45-in.)  diam  tube  with  a  0.25-nim  (0.01-in.)  thick  wall.  The  smallest  test  sections  are  used 
primarily  to  evaluate  the  previously  discussed  wall  material  effects  where  copper  and  its  alloys  are 
much  more  difficult  to  electrically  heat.  Materials  of  interest  include  Amzirc,  OFHC  copper.  Inconel® 
600,  and  304  stainless  steel.  Material  properties  of  these  materials  are  included  in  Appendix  A.  The 
electrical  resistivity  of  the  materials  from  material  property  reference  manuals  was  verified  at  elevated 
temperature  in  the  AEDC  Precision  Measurement  Equipment  Laboratory.  The  largest  diameter  test 
sections  (fabricated  from  304  stainless  steel)  allow  for  assessment  of  cooling  characteristics  on  test 
sections  with  diameters  approximating  actual  arc  nozzle  throat  diameters.  To  prevent  continuous  op¬ 
eration  of  the  rectifiers  at  the  high  power  conditions  required  by  these  two  test  section  sizes,  the  inter¬ 
mediate  size  test  sections  (fabricated  from  304  stainless  steel)  are  used  for  a  bulk  of  the  HTWL 
parametric  studies.  The  test  sections  are  typically  150  mm  (6  in.)  long,  although  shorter  lengths  can 
be  used  to  assess  length  effects.  The  intermediate  size  test  sections  are  fabricated  from  stock  tubing. 
All  other  test  section  tubes  are  fabricated  using  the  EDM  (electrical  discharge  machining)  process  at 
AEDC.  A  uniform  surface  finish  between  the  test  sections  is  maintained  to  prevent  differing  contri¬ 
butions  of  surface  roughness  on  heat  transfer.  Typical  rms  roughness  (measured  with  a  Taylor-Hob- 
son  Surtronic  3P  Surface  Roughness  Machine)  on  the  outside  surface  of  the  tubes  is  approximately 
0.76  to  0.91  lim  (30  to  36  |J.in.).  Wall  thickness  variations  are  determined  with  a  Zeiss  Coordinate 
Measurement  Machine.  Typical  standard  deviation  for  a  24-point  measurement  on  the  19-mm  (0.75- 
in.)  diam  tube  with  a  0.51-mm  (0.02-in.)  thick  wall  is  0.005  mm  (0.0002  in.).  Each  test  section  is  filled 
with  the  low  thermal  conductivity  epoxy,  Sauereisen®  31,  to  prevent  structural  deformation  due  to  the 
high  water  pressure,  and  to  protect  the  internal  instrumentation. 


^For  an  annulus,  the  hydraulic  diameter  is  defined  as  two  times  the  annular  gap,  or  the  outer  diameter  minus  the  inner 
diameter. 
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The  tubular  test  sections  are  press  fit  onto  the  copper  end  pieces  approximately  7.6  mm  (0.3 
in.).  The  initial  method  of  using  high-temperature  silver  solder  to  attach  the  test  section  tubes  to  the 
copper  end  pieces  was  found  to  be  unsuccessful  when  damage  to  internal  instrumentation  and  leaks 
from  the  joints  were  noted  in  early  testing.  Subsequent  joints  made  with  the  AEDC  electron  beam 
welding  technique  yielded  satisfactory  results.  A  hydrostatic  pressure  test  bottle  is  used  to  leak 
check  the  test  section  welds  prior  to  installation  in  the  HTWL  apparatus.  Figure  9a  shows  a  19-mm 
(0.75-in.)  diam,  stainless-steel  test  section  welded  to  the  electrolytic  tough  pitch  copper  end  pieces, 
and  a  closeup  of  the  electron  beam  weld  is  shown  in  Fig.  9b. 

3.4  INSTRUMENTATION  AND  DATA  ACQUISITION 

Generally,  the  most  important  measurements  to  be  made  in  a  flow  boiling  apparatus  are  the 
coolant  conditions,  the  energy  dissipated  at  the  test  section,  and  the  surface  temperature  of  the  test 
section.  Instmmentation  in  the  HTWL  to  monitor  and  record  coolant  temperature,  pressure,  and 
flow  rate  at  various  points  in  the  flow  circuit  allows  determination  of  the  coolant  conditions.  This 
instrumentation  includes  absolute  temperature  (type  T  thermocouples)  and  pressure  (0-2000  psi 
Viatran  pressure  transducers)  measurements  in  the  inlet  and  outlet  coolant  manifolds  and 
differential  temperature  measurements  between  the  manifolds.  Coolant  flow  rate  is  determined 
from  a  51-mm  (2-in.)  diam  Hoffer  Flow  Controls  turbine  flowmeter  installed  upstream  of  the  cool¬ 
ant  pressure  regulator  and  inlet  manifold.  In  addition,  the  total  flow  rate  and  inlet/outlet  pressures 
of  the  high-pressure  demineralized  water  pump  are  monitored  in  the  HTL  main  control  room  during 
an  HTWL  run.  Absolute  pressure  of  the  coolant  is  measured  with  0-2000  psi  Teledyne  Tabor 
pressure  transducers  at  various  axial  locations  (at  the  outside  wall  of  the  annulus)  along  the  length 
of  the  test  section.  Coolant  pressure  drop  (Statham  Pressure  Transducer)  and  a  high-speed  pressure 
measurement  (Kulite®  Pressure  Transducer)  are  also  recorded  at  the  test  section. 

Several  methods  are  used  to  determine  the  power  or  energy  existing  at  the  test  section.  A 
discussion  of  the  actual  approaches  used  in  determining  various  energy  balances  is  included  in 
Section  3.5,  Data  Reduction.  Instmmentation  used  to  support  the  energy  balance  computations 
include  d-c  current  (internal  rectifier  shunt)  and  voltage  for  each  rectifier,  the  total  power  produced 
by  each  rectifier  (measured  with  an  Ohio  Semitronics  Model  PCS  Watt  Transducer),  the  tme  rms 
current  and  voltage  for  each  rectifier  (measured  with  an  Ohio  Semitronics  Model  VT8  Variable 
Frequency  Voltage  Transducer),  and  the  test  section  d-c  voltage  drop.  The  test  section  rms  voltage 
drop  is  proportional  to  the  average  rms  voltage  for  the  rectifiers. 

Probably  the  most  difficult  measurement  is  the  surface  temperature  of  the  test  section.  Ideally, 
a  surface  temperature  at  the  heated  wall  and  coolant  interface  is  desired  for  heat-transfer  analyses. 
In  reality,  an  intmsive  measurement  would  affect  the  wall  heat  transfer  or  dismpt  the  coolant  flow 
pattern,  and  a  non-intmsive  measurement  of  surface  temperature  at  the  interface  is  very  difficult,  if 
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not  impossible,  to  make.  In  the  HTWL  test  sections,  No.  30  (0.25-mm  or  0.01-in.  diam)  type  K 
thermocouples  are  attached  to  the  inside  wall  of  the  test  section  at  various  stations  along  its  length, 
permitting  the  coolant/wall  interface  temperature  (on  the  opposite  side  of  the  wall)  to  be  determined 
analytically.  A  discussion  of  the  inferred  coolant/wall  interface  temperature  computations  is 
included  in  the  next  section  and  Appendix  C.  The  thermocouples  are  typically  spot  welded  to  the 
wall  of  the  larger  diameter  test  sections  and  glued  with  the  high  thermal  conductivity  epoxy,  East¬ 
man  P-10,  in  the  smaller,  thin-walled  (0.25-mm  thick  wall)  test  sections.  In  both  installations,  one 
leg  of  the  thermocouple  is  attached  directly  to  the  surface,  and  the  junction  is  made  approximately 
0.75  mm  (0.03  in.)  above  the  surface,  as  recommended  by  Hughes  in  Ref.  21.  This  arrangement 
prevents  erroneous  temperature  indications  from  a  voltage  drop  produced  by  the  current  flow  in  the 
test  section.  No  separation  of  the  thermocouple  wires  from  the  tube  surface,  caused  by  thermal 
expansion  of  the  test  section,  has  been  noted  in  posttest  inspections.  Steady-state  temperature 
response  from  epoxied  and  spot-welded  thermocouples  on  a  representative  test  section  compared 
within  1.1  °C  (2°F)  in  a  laboratory  oven  up  to  480°C  (900°F).  In  addition,  only  a  slight  conduction 
effect,  due  to  the  presence  of  the  thermocouples,  was  verified  using  the  2-D  axisymmetric  heat 
conduction  program  TRAX  (Ref.  22).  Teflon  insulated  thermocouples  were  used  during  initial  test¬ 
ing;  however,  internal  test  section  temperatures  exceeded  the  vaporization  temperature  of  the 
Teflon,  and  shorting  of  the  thermocouple  wires  was  experienced.  The  problem  was  alleviated  by 
switching  to  braided  glass  insulated  thermocouples,  although,  care  had  to  be  exercised  when  using 
the  thermocouple  wire  in  a  damp  environment.  As  mentioned  previously,  each  test  section  is  filled 
with  the  low  thermal  conductivity  epoxy,  Sauereisen  31,  to  provide  additional  protection  to  the 
internal  thermocouple  wires. 

The  critical  heat  flux  in  a  flow  boiling  apparatus  typically  occurs  at  the  most  downstream 
location  of  the  heated  test  piece,  and  surface  temperature,  therefore,  is  of  primary  interest  at  that 
location.  Conduction  heat-transfer  effects  in  the  tubular  test  sections,  caused  by  the  presence  of  the 
copper  end  pieces,  necessitate  analytical  modeling  of  the  configuration  for  the  determination  of  the 
best  placement  of  the  thermocouples  near  the  ends  of  the  test  section  for  accurate  temperature 
measurement.  The  3-D  thermostructural  computer  code,  ANSYS®,  is  used  to  model  the  conduction 
heat-transfer  effects  in  the  test  sections  of  different  size  and  material  at  various  power  levels.  The 
code  also  allows  for  the  inclusion  of  temperature  sensitive  material  properties  such  as  the  thermal 
and  electrical  conductivities.  Figure  10  shows  a  typical  temperature  distribution  of  an  HTWL  test 
section  with  the  copper  end  pieces.  It  was  found  that  placement  of  the  thermocouples  at  least  10  mm 
(0.4  in.)  from  the  ends  of  the  heated  portion  of  the  test  sections  reduced  the  conduction  effect  to  an 
acceptable  level. 

Additional  instrumentation  in  the  HTWL  includes  type  T  thermocouples  on  each  ballast 
element.  The  thermocouples  are  electrically  isolated  from  the  current  carrying  ballast  elements  with 
thin  mica  sheets.  Type  T  thermocouple  probes  in  the  raw  water  supply  and  return  manifolds  for  the 
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ballast  resistor  bank  provide  cooling  water  temperature  monitoring.  Analog  pressure  gages  located 
in  the  demineralized  water  and  raw  water  manifolds  allow  for  quick  assessment  of  water  pressure 
conditions.  Because  of  the  high  water  pressure  and  high  electrical  energy  present  in  the  apparatus, 
a  modular  control  room  was  installed  on  the  HTL  mezzanine  near  the  HTWL  apparatus  (Figs.  5  and 
11)  to  protect  the  HTWL  data  acquisition  system,  ground  fault  protection  system,  and  personnel. 
Remote  control  of  the  rectifiers  and  operation  of  the  HTWL  flow  control  devices  are  possible  from 
the  control  room.  The  thermocouple  harnesses  are  twisted  (approximately  three  turns  per  inch)  and 
shielded  between  the  apparatus  and  the  control  room  to  reduce  interference  caused  by  the  proximity 
of  the  high  current  flow  in  the  test  section  and  electrical  feed  cables.  In  order  to  further  isolate  the 
test  section  thermocouple  measurements,  each  signal  is  fed  through  a  Preston  Amplifier  prior  to 
sampling  by  the  data  acquisition  system.  Monitoring  of  the  rectifier  performance  (a-c  ripple 
effects)  is  accomplished  with  a  rack-mounted  Tektronix  8300  XWB  Oscilloscope  (Fig.  11)  and  a 
portable  Tektronix  2445  Trigger  Oscilloscope. 

Approximately  60  channels  of  data  are  recorded  during  each  run  with  a  NefP  Instrument 
Corp.  Model  470  data  acquisition  system  (Fig.  11).  Steady-state  data  are  typically  recorded  at  20 
samples  per  sec  for  a  5-sec  burst  at  each  power  setting.  Thermocouple  ice  point  references  and 
system  calibration  are  provided  by  the  Neff,  although  external  ice  point  references  are  required  by 
the  test  section  thermocouples  since  the  Preston  amplifiers  are  used.  The  Neff 470  also  permits  real¬ 
time  monitoring  of  pertinent  measurements  during  a  run.  In  addition,  the  Neff  provides  contact 
closures  based  on  adjustable  parameter  limits  for  use  as  system  interlocks.  Currently,  the  Neff 
initiates  rectifier  shutdown  if  a  low  limit  on  demineralized  water  flow  rate  or  a  high  limit  on  ballast 
resistor  element  temperatures  is  reached.  A  schematic  of  the  data  acquisition  and  monitoring 
systems  is  presented  in  Fig.  12. 

3.5  DATA  REDUCTION 

Raw  data  from  the  Neff  data  acquisition  system  are  averaged  for  each  power  setting  and 
used  as  the  primary  input  for  the  FORTRAN  data  reduction  program.  The  data  reduction  program 
is  used  to  compute  various  energy  balances,  demineralized  water  coolant  conditions,  total  heat 
flux  generated  at  the  HTWL  test  section,  and  the  steady-state  internal  temperature  distribution  in 
the  test  section. 

Several  methods  are  used  to  determine  the  power  or  energy  introduced  to  the  test  section.  One 
method  involves  an  energy  balance  using  the  coolant  mass  flow  and  temperature  rise  measurements 
at  the  test  section.  Another  method  makes  use  of  the  current  flow  and  the  voltage  drop  across  the 
test  section.  A  slight  variation  of  this  method  makes  use  of  the  test  section  material  resistivity  rather 
than  the  test  section  voltage  drop.  The  actual  power  produced  by  each  rectifier  is  measured  with  an 
Ohio  Semitronics  Model  PCS  Watt  Transducer  for  comparison  with  each  of  the  energy  computation 
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methods  mentioned.  The  heat  flux  based  on  the  energy  dissipated  is  then  computed  using  the  test 
section  geometry. 

An  accurate  calculation  of  the  temperature  distribution  in  the  HTWL  test  section  is  required 
to  determine  the  surface  temperature  at  the  test  section  wall/coolant  interface.  Initial  methods  used 
to  compute  the  steady-state  temperature  distribution  included  a  simple  integration  of  the  steady- 
state  1-D  planar  and  radial  conduction  heat-transfer  equations  with  uniform  volumetric  heat 
generation.  Because  the  methods  assume  constant  material  properties  and  little  effect  of  the  extreme 
convection  boundary  condition,  a  more  appropriate  1-D  axisymmetric  finite  difference  approach 
was  chosen  to  achieve  higher  accuracy.  The  approach  allows  for  temperature-dependent  material 
properties  (thermal  and  electrical  conductivity),  and  hence,  nonuniform  internal  heat  generation.  In 
addition,  the  flow  of  current  is  allowed  to  redistribute  within  the  wall  thickness,  depending  on  the 
local  material  resistance.  An  adiabatic  wall  is  assumed  on  the  inside  surface  of  the  test  section 
where  the  low  thermal  conductance  epoxy  is  present  in  the  actual  configuration.  An  initial  estimate 
of  the  convective  heat-transfer  coefficient  on  the  outside  surface  of  the  test  section  is  obtained  from 
boiling  heat-transfer  correlations  such  as  those  discussed  in  Ref.  1 1 .  Gauss-Seidel  iteration  with 
relaxation  is  used  to  reach  a  steady-state  temperature  solution.  Once  a  temperature  distribution 
solution  meeting  the  selected  error  criteria  is  reached,  the  convective  heat-transfer  coefficient  is 
adjusted  until  the  computed  inside  surface  temperature  matches  the  experimentally  measured  wall 
temperature  from  the  test  section  thermocouples.  The  coefficient  adjustment  procedure  may  be  by¬ 
passed  if  the  measured  wall  temperature  is  known  to  be  in  error. 

Because  of  possible  voltage  losses  in  the  actual  HTWL  hardware,  the  input  test  section  voltage 
drop  is  adjusted  such  that  the  final  calculated  total  current  matches  the  measured  current  at  a  given 
power  level  from  the  HTWL  experiment.  An  additional  total  heat  flux  value  is  computed  from 
summation  of  the  individual  element  heat  fluxes  in  the  finite  difference  temperature  calculation. 

A  more  complete  discussion  of  the  data  reduction  equations  is  provided  in  Appendix  C. 

4.0  TEST  PROCEDURE 

Installation  and  removal  of  the  HTWL  test  section  are  relatively  straightforward.  The 
electrical  contact  surfaces  of  the  test  section  copper  end  pieces  (including  the  downstream  piston 
assembly)  are  coated  with  Cool- Amp  Silver  Plating  Powder  and  Bumdy  Penetrox  A  Anti-Oxidant 
to  prevent  arcing  between  the  slip  fit  surfaces.  The  downstream  end  piece  cap  is  then  installed. 
Following  the  cleaning  of  the  test  section  surface  with  denatured  alcohol,  the  two  halves  of  the 
water  jacket  are  assembled  around  the  test  section.  The  assembly  is  then  slipped  into  the  test  section 
assembly  shell  (Fig.  8b),  assuring  that  the  alignment  pin  hole  in  the  water  jacket  is  aligned  with  the 
alignment  pin  port  on  the  shell.  Once  aligned,  the  alignment  pin  is  installed  with  the  appropriate  o- 
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ring.  The  stainless-steel  upstream  cap  and  copper  downstream  cap  are  installed  with  the  appropriate 
o-rings.  Concentricity  of  the  test  section  with  the  water  jacket  is  verified  by  measuring  the  distance 
from  selected  instrumentation  port  faces  to  the  test  section  with  a  depth  gage.  Hookup  of  the  test 
section  thermocouples  to  the  permanent  thermocouple  patch  panel  completes  the  installation. 
Removal  of  the  test  section  upon  test  completion  is  accomplished  in  the  reverse  order. 

The  HTWL  demineralized  water  circuit  is  filled  using  a  57  liters  per  min  (15  gpm)  makeup 
pump  (Union  Pump  Co.  Model  TD-50)  and  venting  air  out  of  the  circuit  through  the  vent  valves  on 
top  of  the  HTWL  test  section  and  inlet  and  outlet  manifolds.  Once  the  circuit  is  filled  (water 
discharge  from  vent  lines),  the  vent  valves  are  closed  and  the  circuit  is  pressurized  to  34.5  bar  (500 
psi),  using  the  makeup  pump  for  leak  check  purposes.  The  ballast  resistor  bank  is  adjusted  to  the 
desired  length  based  on  the  voltage  drop  anticipated  at  the  test  section.  Low-pressure  raw  water 
flow  is  established  at  the  ballast  resistor  bank.  The  HTWL  control  room  power,  rectifier  control 
power,  data  acquisition  system,  and  computer  are  turned  on,  and  instrumentation  calibration  is  ini¬ 
tiated  using  the  Neff  data  acquisition  system. 

A  typical  run  sequence  begins  with  the  manual  startup  of  the  rectifier  cooling  fans.  Water 
flow  through  the  apparatus  is  established  using  the  HTL  high-pressure  demineralized  water  pump. 
The  coolant  flow  rate  and  pressure  at  the  test  section  are  then  adjusted  to  the  planned  test  condition. 
The  ground  fault  panel  is  then  reset  prior  to  rectifier  startup.  Electrical  power  to  the  test  section  is 
initiated  by  closing  the  power  feed  circuit  breaker,  starting  each  rectifier  at  a  low  output  level,  and 
slowly  increasing  power  to  the  first  set  point.  Following  stabilization  of  various  measurements, 
approximately  5  sec  of  data  are  acquired  with  the  Neff  data  acquisition  system,  after  which  the 
power  is  slowly  increased  to  the  next  power  set  point.  At  least  two  power  cycles  are  performed  with 
each  test  section  prior  to  burnout  to  assess  aging  effects  or  identify  other  data  hysteresis.  The 
estimated  CHF  (based  on  pretest  predictions  and  data  acquired  previously  at  the  same  test 
conditions)  is  approached  slowly  with  small  increments  in  power  until  the  test  section  fails  at  burn¬ 
out.  Test  section  failure  is  accompanied  by  a  change  in  power  demand  from  the  rectifiers  which,  in 
turn,  phases  back  the  rectifiers  to  a  negligible  power  setting  prior  to  shutdown.  The  use  of  a  burnout 
detector  to  prevent  test  section  destruction  is  not  possible  because  of  the  speed  and  intensity  of  the 
transition  from  nucleate  to  film  boiling  (Ref.  11).  Demineralized  water  conductivity  and  dissolved 
oxygen  content  levels  are  recorded  prior  to  and  after  each  run  with  an  Omega  Engineering  Inc. 
Model  PHH-10  Conductivity/Temperature/PH  Meter  and  a  Cole  Parmer  Model  5946-70  Dissolved 
Oxygen  Meter,  respectively. 

5.0  SYSTEM  PERFORMANCE  AND  DATA  UNCERTAINTY 

Initial  test  results  obtained  after  shakedown  and  checkout  of  the  HTWL  aided  in  the  definition 
of  the  performance  envelope  and  operating  characteristics  of  various  support  systems.  Figure  13 
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presents  the  flow  performance  envelope  for  the  three  test  sections  discussed  in  Section  3.3.  The 
available  volumetric  flow  rate  at  a  given  inlet  manifold  pressure  measured  just  upstream  of  the 
HTWL  test  section  is  shown  in  Fig.  13a  and  at  a  given  test  section  exit  pressure  (burnout  location) 
in  Fig.  13b.  The  water  jackets  selected  for  each  test  section  provided  for  hydraulic  diameters  of  3.76 
mm  (0.148  in.)  for  the  26-mm  (1.024-in.)  diam  tube,  5.08  mm  (0.2  in.)  for  the  19-mm  (0.75-in.) 
diam  tube,  and  7.06  mm  (0.278  in.)  for  the  12-mm  (0.45-in.)  diam  tube.  The  optimum  setting  for 
the  demineralized  water  bypass  valve,  which  resulted  in  the  highest  mass  flow  achievable,  yet 
providing  for  safe  starting  and  operation  of  the  high-pressure  pump,  was  found  to  be  one-fifth  open. 
Generally,  the  highest  output  pressures  from  the  HTL  high-pressure  demineralized  water  pump  are 
achieved  with  a  makeup  pump  suction  pressure  of  48.3  bar  (700  psi).  Higher  pressures  are  attain¬ 
able  at  the  higher  flow  rates;  however,  the  pump  is  limited  to  operation  below  124  bar  (1800  psi)  at 
the  pump  discharge  or  approximately  114  bar  (1650  psi)  at  the  inlet  manifold. 

Rectifier  performance  is  presented  in  Fig.  14.  Each  of  the  two  rectifiers  was  operated  in  the 
current  control  mode  up  to  100-percent  power  output  by  installing  a  solid  copper  test  section  and 
using  the  maximum  available  ballast  resistor  length.  Shown  in  Fig.  14a  are  the  d-c  current  output 
and  rms  current  output  (the  actual  output  of  each  rectifier  made  up  of  the  d-c  signal  with  an  a-c  com¬ 
ponent  riding  on  the  d-c  waveform)  for  each  rectifier.  Also  included  in  the  figure  are  the  combined 
d-c  and  rms  current  outputs  of  the  two  rectifiers  since  power  requirements  typically  necessitated  the 
parallel  operation  of  the  rectifiers.  Figure  14b  shows  the  power  output  of  the  rectifiers  up  to  100 
percent.  Voltage  limitation  of  the  ballast  resistor  bank  prevented  the  output  of  the  full  1 .6  MW  of 
power  at  100  percent  for  the  solid  test  section.  As  mentioned  previously,  the  rectifiers  may  be  op¬ 
erated  up  to  200  percent  for  limited  mn  times. 

The  ballast  resistor  voltage  and  surface  temperature  characteristics  as  a  function  of  rectifier 
rms  current  are  presented  in  Fig.  15.  As  in  the  rectifier  performance  checks,  a  solid  copper  test 
section  with  a  negligible  voltage  drop  was  used  to  obtain  the  ballast  resistor  characteristics.  The 
maximum  voltage  drop  that  can  be  obtained  across  the  ballast  resistor  elements  is  80  volts  (Fig. 
15a).  Surface  temperatures  on  the  elements  remain  below  140°C  (300°F),  as  shown  in  Fig.  15b. 

Electromagnetic  fields  generated  by  the  high  power  electrical  systems  were  of  concern  for 
personnel  safety  and  control  room  equipment  reliability,  and  were  evaluated  during  initial  testing 
in  HTWL.  Magnetic  flux  density  (MFD)  measurements  were  performed  with  an  F.  W.  Bell  Model- 
9500  Gaussmeter  in  order  to  quantify  the  level  of  electromagnetic  interference  in  the  HTWL  control 
room  during  rectifier  operation.  Figure  16  presents  the  a-c  (unfiltered)  and  d-c  MFD  measurements 
as  a  function  of  total  rectifier  power.  The  a-c  MFD  did  not  change  with  rectifier  operation  or 
measurement  location  in  the  control  room;  however,  d-c  MFD  measurements  showed  considerable 
effects  of  rectifier  power  level  and  location  of  the  measurement  within  the  control  room.  The 
maximum  d-c  MFD  levels  were  recorded  near  the  floor  of  the  control  room,  and  are  shown  in  Fig. 
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16  along  with  the  maximum  values  recorded  at  chest  height  throughout  the  room.  No  differences  in 
the  MFD  levels  were  detected  with  the  absence  of  6900-V  power  fed  to  the  HTL  building  and  with 
6900-V  power  to  the  rectifiers  with  zero  output  (i.e.,  no  change  in  MFD  levels  with  or  without  the 
presence  of  the  6900-V  power  feed). 

Figure  17  presents  typical  data  acquired  in  the  HTWL  for  a  19-mm  (0.75-in)  diam,  304 
stainless-steel  test  section  at  the  noted  test  conditions.  The  curve  presented  in  Fig.  17a  is  a  pretest 
prediction  of  the  boiling  curve  using  Kays  and  Leung  (Ref.  23)  correlation  for  pure  forced 
convection,  Bergles  and  Rohsenow  (Ref.  4)  correlation  for  transition  from  pure  convection  to  fully 
developed  nucleate  boiling,  and  the  Rohsenow  (Ref.  5)  nucleate  boiling  correlation  in  the  boiling 
regime.  As  can  be  seen  in  Fig.  17a,  lower  wall  temperatures  have  been  noted  in  the  experiments 
than  are  predicted  by  the  above  correlations;  however,  good  agreement  is  shown  near  burnout.  A 
possible  cause  for  the  discrepancy  may  be  the  inaccuracy  in  the  wall  temperature  measurement.  As 
noted  previously,  care  had  to  be  exercised  in  preventing  the  braided  glass  insulation  from  becoming 
wet.  The  presence  of  a  small  amount  of  moisture  on  the  thermocouple  leads  resulted  in  considerable 
noise  in  the  temperature  measurements.  A  possible  correction  for  the  moisture  problem  would  be 
to  incorporate  metallic-sheathed  thermocouple  leads.  The  CHF  predictions  from  the  correlations  of 
Bemath  (Ref.  12),  Van  Huff  and  Rousar  (Ref.  13),  Rousar  (Ref.  14),  Yagov  and  Puzin  (Ref.  15), 
Levy  (Ref.  16),  and  Labuntsov  (Ref.  17)  are  also  shown  in  Fig.  17a.  The  CHF  occurred  at  the  most 
downstream  station  of  the  test  section,  and  therefore,  the  data  presented  in  the  figure  were  measured 
at  that  particular  station.  The  correlations  of  Labuntsov  and  Rousar  best  predicted  the  CHF  in  this 
particular  test;  however,  neither  correlation  included  data  from  annular  configurations. 

Figure  17b  shows  a  comparison  of  the  various  heat  flux  calculations.  The  total  heat  flux  (qrox) 
from  the  hardware  is  determined  from  the  rms  current  and  voltage  measured  at  the  rectifiers  and, 
therefore,  has  the  largest  value  because  of  energy  losses  in  the  electrical  cables,  ballast  resistor,  and 
attachments.  The  test  section  measured  heat  flux  (qyg)  is  determined  similarly  to  the  total  heat  flux, 
except  that  the  rms  voltage  drop  across  the  test  section  assembly  is  used  instead  of  the  rectifier  voltage. 
The  heat  flux  computed  from  the  coolant  temperature  rise  through  the  test  section  assembly  (q  sys)  and 
the  heat  flux  calculated  from  a  summation  of  the  element  internal  heat  generation  in  the  finite 
difference  routine  (q  calc)  agree  within  6  percent  at  all  but  the  lowest  power  settings. 

Typical  dependance  of  the  pressure  drop  across  the  test  section  assembly  on  heat  flux  is 
illustrated  in  Fig.  17c.  As  heat  flux  is  initially  increased,  the  pressure  drop  decreases  because  of 
decreasing  friction  factor  (see  Ref.  24).  As  boiling  begins  and  becomes  well  established,  the 
pressure  drop  increases.  As  shown  in  Fig.  17c,  the  pressure  drop  across  the  150-mm  (6-in.)  long  test 
section  is  small  at  the  elevated  pressure  and  velocity  for  the  particular  test. 
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The  effects  of  boiling  at  elevated  heat  flux  has  not  been  detected  in  the  form  of  pressure 
oscillations.  Figure  17d  presents  typical  high-speed  pressure  measurements  for  low  power  settings 
where  no  boiling  could  exist  and  higher  settings  where  boiling  was  suspected.  Only  slight 
differences  between  the  nonboiling  and  boiling  results  can  be  seen  in  the  figure. 

Burnout  of  the  test  sections  typically  occurred  at  the  most  downstream  location  of  the  heated 
tube  as  shown  in  Fig.  18.  Because  the  bulk  fluid  temperature  is  highest  and  the  thermal  boundary 
layer  is  largest  at  this  location,  it  follows  that  the  cooling  would  be  less  efficient,  thereby  promoting 
burnout.  However,  during  preliminary  testing  in  the  HTWL,  a  few  burnouts  occurred  at  the  most 
upstream  location  on  the  test  pieces.  Typically,  the  upstream  burnouts  occurred  at  lower  power 
settings  indicating  premature  failure  due  to  structural  anomalies  or  an  inadequate  weld.  Such 
problems  would  cause  coolant  leakage  to  the  interior  of  the  test  piece,  deterioration  of  the  internal 
epoxy  support,  and  eventual  collapse  and  melting  of  the  heated  tube  wall. 

The  rectifier  ripple,  or  the  a-c  component  that  rides  on  the  d-c  output  signal,  contributes  a 
significant  amount  to  the  uncertainty  of  the  instantaneous  power  output  of  the  rectifiers.  The  ripple 
for  the  two  rectifiers  was  measured  during  the  rectifier  performance  checks  using  the  solid  test 
section  described  previously.  The  measured  ripple  for  each  rectifier  is  shown  in  Fig.  19a,  and  the 
actual  waveforms  are  presented  in  Fig.  19b.  Because  the  rectifiers  are  operated  near  the  maximum 
voltage  output  of  100  volts,  the  ripple  accounts  for  approximately  5  to  7  percent  uncertainty  in  the 
instantaneous  power  output.  Use  of  rms  values  for  rectifier  current  and  voltage  in  the  dissipated 
energy  and  heat  flux  equations  minimizes  the  effects  of  ripple  on  the  time-averaged  data 
uncertainty.  Uncertainty  in  the  rms  current  and  voltage  is  primarily  measurement  device  and  data 
system  inaccuracies. 

Additional  factors  which  contribute  to  the  uncertainty  in  the  energy  dissipated  at  the  test 
section  include  variations  in  the  test  section  wall  thickness  and  material  properties.  The  wall 
thickness  variations  are  determined  by  the  AEDC  Precision  Inspection  Laboratory,  as  described  in 
Section  3.3.  The  test  section  material  density  and  electrical  resistivity  from  literature  were  verified 
by  personnel  in  the  AEDC  Precision  Measurement  Equipment  Laboratory.  Uncertainty  of  the  test 
section  wall  thickness  and  material  properties  along  with  measurement  uncertainty  of  other 
instrumentation  used  in  the  HTWL  are  included  in  Table  la. 

In  general,  instrumentation  calibration  and  data  uncertainty  estimates  were  made  using 
methods  recognized  by  the  National  Institute  of  Standards  and  Technology  (NIST).  Measurement 
uncertainty  is  a  combination  of  bias  and  precision  errors  defined  as: 

U,„  =  ±[(B)=  +  (2S)Vn]'“ 


19 


AEDC-TR-95-15 


where  Urss  is  the  root  sum  of  squares  uncertainty,  B  is  the  bias  limit,  S  is  the  standard  deviation 
about  a  mean  value  of  the  measurement  process,  “n”  is  the  number  of  experiment  periods  from 
which  the  samples  were  used  in  determining  the  mean  value  (taken  here  to  equal  one),  and  the 
multiplier  “2”  assumes  10  or  more  samples  associated  with  S  and  is  used  to  ensure  a  95-percent 
coverage  for  the  uncertainty  limits. 

In  addition  to  the  uncertainty,  the  type  and  range  of  measuring  device,  the  type  of  recording 
device,  and  the  method  of  calibration  for  each  measurement  are  provided  in  Table  1.  Propagation 
of  the  bias  and  precision  errors  of  measured  data  through  the  calculated  data  was  made  in 
accordance  with  Ref.  25,  and  the  results  are  presented  in  Table  2. 

6.0  SUMMARY 

In  conclusion,  a  new  high  heat  flux,  flow  boiling  apparatus  for  the  study  of  cooling 
effectiveness  has  been  developed  at  the  AEDC.  The  current  application  of  HTWL  is  for  the 
evaluation  of  cooling  processes  encountered  in  high-pressure,  electric  arc  heater  facilities, 
particularly  arc  heater  nozzles.  The  facility  is  capable  of  providing  up  to  3.2  MW  of  power  to  a 
metallic  test  section,  simulating  heat  flux  levels  in  excess  of  2  x  10*  W/m^  (18,000  Btu/fP  sec), 
with  a  water  flow  rate  up  to  9  kg/sec  (20  Ibm/sec)  at  pressure  up  to  100  bar  (1500  psi).  The  current 
test  section  configuration  allows  for  the  determination  of  parametric  effects  in  annular  flow. 
Early  heat-transfer  data  from  HTWL  show  some  disagreement  with  pretest  predictions  in  the 
fully  developed  nucleate  boiling  regime  and  at  CHF.  Additional  testing  in  the  HTWL  has  been 
performed,  and  reduced  data  are  available  under  separate  cover  from  AEDC.  Future  analyses  of 
the  data  are  planned  so  that  more  accurate  prediction  of  cooling  requirements  for  actual  arc  heater 
nozzles  may  be  accomplished. 
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Water  Cooling  Flow 


High  Enthalpy  Air  Flow 


Figure  1.  Arc  heater  nozzle  heating/cooling. 
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Figure  2.  HTWL  system  schematic. 
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a.  Outside  diameter  of  1  in.  (25.4  mm) 

Figure  3.  Power  supply  requirements  for  an  electrically  heated  stainless-steel  tube. 
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Wall  Thickness,  mm 


b.  Outside  diameter  of  0.75  in.  (19.1  mm) 
Figure  3.  Continued. 
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c.  Outside  diameter  of  0.5  in.  (: 
Figure  3.  Concluded. 
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Figure  4.  Continued. 
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Figure  7.  HTL  high-pressure  demineralized  water  pump. 
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a.  HTWL  test  section 
Figure  8.  HTWL  test  section  assembly. 
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a.  19-mm-diain  stainless-steel  test  section  and  water  jacket 
Figure  9.  HTWL  test  section. 
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b.  Closeup  of  electron  beam  weld 
Figure  9.  Concluded. 
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Figure  10.  Theoretical  temperature  distribution  of  HTWL  test  section. 


AEDC-TR-95-15 


UJ 


Ground  Fault 
Panel 


EUD 

Monitor 


Demineralized  Water 
Control  Panel 


NefT  Data 
Acquisitior 
System 


Rectifier 
Control  Panel 


Instrumentation  Power 
Supplies  and  Amplifiers 


Figure  11.  HTWL  control  room. 
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Figure  12.  HTWL  data  acquisition  and  monitoring  systems. 
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a.  Inlet  manifold  conditions 

Figure  13.  HTWL  demineralized  water  system  performance. 
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b.  Test  section  exit  conditions 
Figure  13.  Concluded. 
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a.  Current  characteristics 
Figure  14.  Rectifier  output  characteristics. 
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b.  Power  based  on  Irms 
Figure  14.  Concluded. 
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Figure  15.  Ballast  resistor  characteristics. 
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b.  Highest  surface  temperature 
Figure  15.  Concluded. 
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Figure  16.  Magnetic  flux  density  measurements  in  HTWL  control  room. 
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a.  Boiling  curve  at  burnout  location  (test  section  exit) 
Figure  17.  Typical  HTWL  test  results. 
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b.  Comparison  of  heat  flux  computation  methods 
Figure  17.  Continued. 
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c.  Pressure  drop  across  test  section 
Figure  17.  Continued. 
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d.  Test  section  high-speed  pressure  response 
Figure  17.  Concluded. 
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Figure  18.  Posttest  HTWL  test  section. 
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a.  Ripple  as  a  function  of  rms  current 
Figure  19.  Rectifier  ripple  characteristics. 
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b.  Rectifier  ripple  waveforms 
Figure  19.  Continued. 
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b.  Concluded 
Figure  19.  Concluded. 
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Table  1.  Basic  Measurement  Uncertainty 


Basic  Measurement 

Range 

Precision 
Limit  2S* 

Bias 
Limit  B* 

Uncertainty 

Urss* 

Type  of  Measuring 
Device 

Type  of 
Recording 
Device 

Method  of  Calibration 

Manifold  Pressure 

13.8-138  bar 
(200  -  2000  psig) 

±0.06 

±0.13 

±0.14 

Viatran 

Transducer 

Neff  Digital 
Data 

Acquisition 

System 

In-place  application  of 
multiple  pressure  levels 
measured  with  a  pressure 
measuring 

device  calibrated  in  the 
Precision  Measure-ment 
Equipment  Lab 

Test  Section 
Pressure 

13.8-138  bar 
(200  -  2000  psig) 

±0.04 

±0.17 

±0.18 

Teledyne  Tabor 
Transducer 

w 

Test  Section 
Pressure  Drop 

1  -  34.5  bar 
(15-500  psi) 

±0.7 

±0.37 

±0.79 

Statham 

Transducer 

u 

Test  Section  High' 
Speed  Pressure 

13.8-138  bar 
(200  -  2000  psia) 

±0.18 

±0.1 

±0.21 

Kulite'^T  ransducer 

u 

Manifold 

Temperature 

10-120°C 
(50  -  250°C) 

±2.4 

±0.72 

±2.5 

Type  T 

Thermocouple 

Neff  Digital 
Data 

Acquisition 

Sytem 

with  Internal 
Reference 

Thermocouple  Verification 
of  NIST  Conformity 

Voltage  Substitution 
Calibration 

Manifold 

Temperature  Drop 

0-110°C 
(32  -  230“C) 

±1.15 

±0.35 

±1.2 

u 

Test  Section 
Temperature 

10-1400”C 
(50  -  2550''F) 

±3.07 

±1.3 

±3.33 

Type  K 

Thermocouples 

NeffDDAS 

with 

External  tee 
Point  Ref.  and 
Preston 
Amp/Filter 

M 

Coolant  Flow  Rate 

0-850  Ipm 
(0  -  225  gpm) 

±1.3 

±0.32 

±1.34 

Hoffer  Flowmeter 

NeffDDAS 

NIST  Conformity  Rate 
Verification  in  the 

Precision  Measurement 
Equipment  Lab 

Current-d'C 

±0.62 

±1.9 

Rectifier  Shunt 

Manufacturer  Calibration: 
In-place  mv  Substitution 
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Table  1.  Concluded 


.  Current-rms 

0  - 10,000  amps 

±1.0 

±2.25 

±2.46 

Ohio  Semitronics 
Voltage 
Transducer 

Neff  DBAS 

u 

Voltage  -  d-c 

0-100  volts 

±1.65 

±0.71 

±1.8 

Rectifier  Bus 

u 

Voltage-rms 

0-100  volts 

±2 

±0.71 

±2.1 

Ohio  Semitronics 
Voltage 
Transducer 

M 

Test  Section  Voltage 
Drop  -  d-c 

0-100  volts 

±1.65 

±0.71 

±1.8 

Volt  Meter 

Rectifier  Pow/er-rms 

S-1600kW 

±2 

0* 

±2 

Ohio  Semitronics 
Watt  Transducer 

Digital  Display 
Only 

Magnetic  Flux 
Density 

d-c;  30  mG  to  30G 
a-c:  50  to  99  Hz 

±0.15 

±3.4 

0* 

±0.15 

±3.4 

F.  W.  Bell 
Gaussmeter 

4( 

Unit  Installed  Standard 

Test  Section 
Roughness 

0  - 10  pm 
(0  -  200  pin.) 

m 

0* 

urn 

Taylor-Hobson 

Roughness 

Machine 

Manufacturer  Calibration 

i 

Test  Section 
Thickness 

0.25  - 1  mm 
(0.01  -  0.04  in.)* 

m 

0* 

±1.0 

Zeiss  Coordinate 
Measurement 
Machine 

M 

^mm 

1 

±1.0 

0* 

±0.1 

Dial  Calipers 

Manual 

Material  Electrical 
Resistivity 

0-140  pohm-cm 
(0  -  -4.5  pohm-ft) 

±2.6 

0* 

±2.6 

Rubicon  Standard 
Resistor/  Voltage 
Measurement 

u 

Comparison  to  Standard 
Resistor  and  Reference 
Data 

*Numerical  values  are  in  “percent  of  reading.' 
+Assumed  to  be  zero 


Table  2.  Calculated  Parameter  Uncertainty 


Parameter 

Precision  Limit,  2S 
%  of  Calculated  Value 

Bias  Limit,  B 
%  of  Calculated  Value 

Uncertainty,  Urss 
%  of  Calculated  Value 

<iTOT»  Etot 

+2 

±0.71 

±2.1 

^SYS»  ^SYS 

±1 .4 

±0.5 

±1.49 

^CALC»  ^CALC 

±3.2 

±0.71 

±3.3 

^TS»  ^TS 

±2 

±0.71 

±2.1 

Coolant  Velocity,  Mass  Velocity 

±1.3 

±0.32 

±1.34 

Calculated  Wall  Temperature 

±3.2 

±0.71 

±3.3 

Test  Section  Voltage  Drop-rms 

±2 

±0.71 

±2.1 
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TEST  SECTION  MATERIAL  PROPERTIES 
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Figure  A-1.  Test  section  material  properties. 
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Figure  A-1.  Concluded. 
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APPENDIX  B 

POWER  COMPUTATION  PROCEDURE 

Compute  the  surface  area  Ag  and  the  cross-sectional  area  Acs  based  on  selected  tube 
diameter,  length,  and  wall  thickness: 

A  ,  in 

r,  da,  L  in  m 
Ags  in 
d^„  di  in  m 

Compute  the  energy  required  to  reach  the  desired  heat  flux  from: 


As  “  Q/TZkIj  —  TzdgL 


E  =  q  A, 


q  inkW/m^ 
As  in 
EinkW 


Compute  the  resistance  of  the  tubular  test  section: 


Pe  in  ohm-m 
L  in  m 

Acs  in 

R  in  ohm 


Compute  the  current  from: 


1 


R  in  ohm 
1  in  amps 

Compute  voltage  from: 

y  =-j  E  in  W 

I  in  amps 
V  in  volts 
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APPENDIX  C 

DATA  REDUCTION  COMPUTATIONAL  APPROACH 

The  primary  calculations  performed  in  the  data  reduction  program  include  the  energy 
balance/heat  flux  computations  and  the  radial  temperature  distribution  in  the  heated  test  section 
(for  the  wall/coolant  interface  temperature  determination). 

Energy/Heat  Flux 

The  total  rectifier  energy  is  calculated  from: 

^TOT  ~  Ams.TOT  ^rms,  avg 

where  Ims  TOT  is  the  sum  of  the  rms  current  output  for  the  two  rectifiers,  and  VREc,avg  is  the  average 
voltage  output  of  the  rectifiers. 

The  corresponding  heat  flux  is  then: 

_  ^TOT 
^TOT  -  ^ 

where  As  is  the  outside  surface  area  of  the  test  section. 

A  somewhat  more  accurate  test  section  result  may  be  obtained  by  using  the  measured 
voltage  drop  across  the  test  section  assembly  rather  than  the  rectifier  voltage,  thereby  avoiding 
power  cable,  ballast  resistor,  and  attachment  losses: 

^TS  ~  4ms,TOT  ^TS.rms 

The  corresponding  heat  flux  is: 
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The  heat  flux  at  the  test  section  may  also  be  determined  by  performing  an  energy  balance  on 
the  coolant  as  it  proceeds  through  the  test  section  assembly.  For  an  open  system  with  no  work  and 
uniform,  steady-state  flow  with  negligible  changes  in  kinetic  and  potential  energy,  the  conserva¬ 
tion  of  energy  equation  reduces  to: 


^SYS  -  ^  (^exil  ~  ^inlet) 

where  m  is  the  mass  flow  of  the  coolant,  and  the  enthalpy  change  for  an  incompressible  substance 
is  given  as: 


^exit  ^inlel 


=  C 


avg 


(  T  exit  -  T inlet  )  +  p  (Pexil  “  ^ inlet  ) 


where  c^v^  is  the  average  specific  heat  of  the  coolant,  p  is  the  coolant  density,  and  T  and  P 
are  the  temperature  and  pressure,  respectively,  at  the  inlet  or  exit  of  the  test  section.  The  contribu¬ 
tion  of  the  pressure  drop  on  the  enthalpy  change  is  small  compared  to  that  of  the  temperature 
change;  therefore: 


and 


^exit  ^inlet 


“  Cavg^T^exit-Ti 


inlet 


-SYS 


-  me 


avg 


i^exii  '^inlet^ 


The  heat  flux  from  the  test  section  then  follows  from: 


%YS  ~ 


^SYS 


The  fourth  approach  for  determining  heat  flux  at  the  test  section  requires  the  use  of  the  test 
section  voltage  drop  and  the  material  electrical  resistivity.  Because  the  resistivity  varies  signifi¬ 
cantly  with  local  temperature  and  the  presence  of  a  severe  radial  temperature  gradient  in  the  test 
section  during  heating,  the  heat  flux  based  on  electrical  resistivity  is  determined  as  part  of  the  test 
section  temperature  distribution  computation  presented  in  the  following  discussion. 


71 


AEDC-TR-95-15 


Radial  Temperature  Distribution 

Because  of  the  severe  radial  temperature  gradient  in  the  test  section,  the  electrical  resistivity 
will  also  vary  in  the  radial  direction,  and  the  internal  heat  generation  will  be  radially  nonuniform 
on  a  per  volume  basis.  Therefore,  a  reasonable  temperature  distribution  is  not  only  important  in 
estimating  the  outside  wall  temperature,  but  is  necessary  in  determining  the  internal  heat 
generation  which,  in  turn,  affects  the  temperature  distribution. 

Using  the  energy  balance  method  and  assuming  no  transverse  heat  conduction,  a  finite 
difference  equation  for  a  node  at  the  center  of  an  internal  element  of  the  test  section  tube  may  be 
obtained: 


q,-_  1  As,  outer  +  1 A5  —  0 


(C-1) 


where 

q,.;  is  the  heat  flux  into  (or  out  of)  the  element  at  the  outer  face, 

^s.outer  is  the  surface  area  of  the  outer  face, 

q,+/  is  the  heat  flux  into  (or  out  of)  the  element  at  the  inner  face, 

^s, inner  is  the  surfacc  area  of  the  inner  face, 

q^  is  the  internal  heat  generation,  and 

Vol  is  the  element  volume. 

The  cylindrical  geometry  of  the  test  section  yields: 

^S, outer  =  (2Ri  +  dr) 

^S.inner  ~  ”  dr) 

Acs  ^  27x7? ,■  dr 
Vol  =  2nR;  drL 
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where  L  is  the  length  of  the  test  section  tube, 

Rj  is  the  radial  distance  to  the  node  i, 
dr  is  the  radial  thickness  of  the  element,  and 
Acs  is  element  cross-sectional  area. 

Invoking  Fourier's  law  of  heat  conduction,  the  heat  fluxes  can  be  expressed  as: 


where  ki+,,  kj.,,  and  k;  are  the  material  thermal  conductivities  at  the  inner  face,  outer  face,  and 
interior  node,  respectively.  Similarly,  T,+i,  Tj.,,  and  Tj  are  the  respective  temperatures.  Substituting 
these  equations  into  Eq.  (C-1)  and  rearranging,  an  equation  for  the  temperature  at  the  interior 
node,  Tj,  is  obtained: 


The  energy  balance  for  the  inner  wall  element  is  treated  in  a  similar  fashion  but  allowing  for 
an  adiabatic  wall  at  the  inner  face: 


where 


l'^^,  outer 


^1+  1^5. 


+  qyol  =  0 


=0 


(C-3) 


^s.outer  ~  (2^/  +  dr), 

A„  =  KRjdr, 

Vol  =  7tRj(irL, 
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and  Ri  is  the  test  section  inner  diameter. 
From  Fourier’s  law: 


Substituting  into  Eq.  (C-3)  and  solving  for  the  inner  wall  temperature,  T,: 

r,  =  [(*,- ,  +  +  i)]  [(*-1  +  (»+  1  +  «/'■'']]  <c-4) 

The  energy  balance  on  the  outer  wall  element  yields: 

l^S,  outer  +  +  1^5.  inner  +  =  0  (C-5) 

where 

^S, outer  ~ 

As.imer  =  TtL(2R,,  -  dr), 

Acj  =  nR^dr, 

Vol  =  nR^drL, 

and  Ro  is  the  test  section  outer  diameter.  From  Fourier's  law  the  inner  face  heat  flux  is: 


where  T,,  is  the  outer  wall  temperature.  From  Newton's  law  of  cooling  the  outer  wall  heat  flux  is 
given  as: 

=  h{Tf^-T„) 

where  h  is  the  assumed  convective  heat-transfer  coefficient,  and  Tf,  is  the  bulk  coolant  temperature. 
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Substituting  into  Eq.  (C-5)  and  solving  for  the  outer  wall  temperature,  T;,: 


'4/?„  JJ 


(C-6) 


The  heat  generation  qg  in  Eqs.  (C-2),  (C-4),  and  (C-6)  is  given  as: 

E 

“  Vol 

However,  the  energy  E  can  be  written  in  the  form  of  electrical  measurements  and  geometry: 


where  E  is  the  generated  energy, 

I  is  the  current, 

V  is  the  voltage  drop  across  the  test  section. 


Therefore, 


Pel 

R  is  the  electrical  resistance  = - ,  and 

Pe  is  the  material  electrical  resistivity. 


p  LVol 


or  simplifying. 


(C-7) 
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The  radial  temperature  distribution  can  be  solved  with  Gauss-Seidel  iteration  by  using  an 
initial  estimate  for  the  temperature  distribution  and  heat-transfer  coefficient,  and  solving  Eqs. 
(C-2),  (C-4),  and  (C-6)  for  a  refined  temperature  distribution  solution.  Substitution  of  this 
refined  solution  into  the  equations  continues  until  a  stable  solution  is  achieved  (a  small  residual 
between  successive  temperature  calculations).  The  number  of  iterative  cycles  can  be  reduced  by 
using  relaxation  of  the  form: 


T  =37  ,+  (1  -  3)  r 

new  ^  current  ^  i 


previous 


where  Tcu„cnt  is  the  temperature  computed  from  Eqs.  (C-2),  (C-4),  and  (C-6); 

Tprevious  is  the  temperature  previously  calculated  and  used  in  the  equations 
to  compute  T,„„ent 

3  is  the  relaxation  factor  ;  for  linear  systems  1  <  3  <  2; 

Tnew  is  the  new  temperature  to  be  used  in  the  next  iterative  cycle. 


The  heat  rates  from  each  elemental  heat  generation  may  be  summed  for  all  elements  when 
temperature  convergence  is  achieved,  and  the  total  heat  flux  from  the  test  section,  4  calc  is 
calculated  using  the  overall  test  section  geometry.  The  corresponding  calculated  energy  is 

^CALC  ~  ^CALC  A 

Since  the  inner  wall  temperature  is  measured  at  several  stations  along  the  test  section  length, 
a  direct  comparison  between  the  measured  value  and  the  calculated  inner  wall  temperature  can  be 
made.  Some  difference  between  the  two  temperatures  is  expected  because  the  convective  heat- 
transfer  coefficient  on  the  outer  surface  of  the  test  section  that  is  used  in  the  computations  is  only 
an  estimate.  This  heat-transfer  coefficient  can  be  adjusted  within  the  data  reduction  program  to 
provide  a  match  between  the  measured  and  calculated  values  of  the  inner  wall  temperature. 

In  addition  to  the  heat  rates  from  each  elemental  heat  generation,  the  current  flow  through 
each  element  may  be  summed  to  obtain  the  total  current  flow  through  the  test  section.  Because  the 
measured  voltage  drop  across  the  test  section  has  some  uncertainty  associated  with  a-c  ripple  and 
losses  in  slip  fit  joints  within  the  test  section  assembly,  the  voltage  drop  used  in  the  heat  generation 
calculation,  Eq.  (C-7),  may  be  adjusted  within  the  program  such  that  the  calculated  current  and 
measured  rms  current  match.  Such  voltage  drop  corrections  are  typically  small. 


The  selection  of  the  number  of  elements  in  the  finite  difference  scheme  is  important  such 
that  continuity  is  maintained;  i.e.,  the  resistance  to  conduction  within  the  test  section  is  much  less 
than  the  resistance  to  convection  across  the  coolant  boundary  layer.  The  thickness,  dr,  of  the 
elements  is  selected  such  that  the  Biot  number,  Bi,  remains  less  than  0.1: 
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<0.1 


where  h  is  the  largest  heat-transfer  coefficient  expected,  and  k  is  the  lowest  test  section  material 
thermal  conductivity  expected. 

For  comparison  purposes,  a  simplistic  solution  of  the  integrated  cylindrical  1-D  heat-transfer 
equation  can  be  used  to  compute  a  test  section  outer  wall  temperature.  Assuming  steady-state 
conditions,  constant  material  properties,  uniform  heat  generation,  and  no  transverse  heat 
conduction,  the  heat  equation  reduces  to 


rdr\  dr) 


=  0. 


For  the  HTWL  test  section  the  appropriate  boundary  conditions  are: 

1 .  T  (Ri)  =  Ti  measured  with  a  thermocouple 
dT  I 

2.  -—  /?,  =  0  assumed  adiabatic  inner  wall 
dr\ 

Integrating  the  above  heat  equation  and  invoking  the  boundary  conditions,  the  temperature  as  a 
function  of  radial  position,  r,  is 


and  solving  for  the  outer  wall  temperature: 

To  =  +  (C-8) 

Because  of  the  severe  temperature  gradients  expected  and,  hence,  the  nonuniform  material 
properties  and  heat  generation  per  unit  volume  in  the  radial  direction,  Eq.  (C-8)  is  expected  to  have 
a  significant  error. 
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NOMENCLATURE 

Acs  Test  section  cross-sectional  area,  or  ft^ 

As  Test  section  outside  surface  area,  or  ft^ 

B  Bias  limit 

CHF  Critical  heat  flux,  W/m^  or  Btu/fF  sec 

Dh  Hydraulic  diameter,  mm  or  in. 

dTsuh  Degree  of  subcooling,  °C  or  °F  (Ts^,  -  Tb) 

E  Energy,  kW 

G  Mass  velocity,  kg/m^  sec  or  Ibm/fF  sec 

h  Heat-transfer  coefficient,  W/m^  °C  or  Btu/fF  sec  °F 

HTL  High  Temperature  Laboratory 

HTWL  High  Temperature  Wall  Laboratory 
I  Current,  amps 

d-c  current,  amps 
Irnis  rms  current,  amps 

L  Test  section  length,  m  or  ft 

Pexit  Test  section  exit  pressure,  bar  or  psi 

Psuction  Back  pressure  on  the  high-pressure  demineralized  water  pump,  bar  or  psi 
q  Heat  flux,  W/m^  or  Btu/fF  sec 
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q  CALC  Calculated  heat  flux  from  temperature  distribution  finite  difference  routine, 
W/m^  or  Btu/ff  sec 

q  SYS  Heat  flux  computed  from  temperature  rise  of  coolant,  W/m^  or  Btu/ff  sec 

q  tot  Total  heat  flux  computed  from  rectifier  current  and  voltage,  W/m^  or  Btu/ff  sec 

q  Ts  Heat  flux  computed  from  rectifier  current  and  test  section  voltage  drop,  W/m^ 

or  Btu/ff  sec 

R  Resistance,  ohms 

S  Standard  deviation 

Tb  Bulk  coolant  temperature,  °C  or  °F 

Tsai  Coolant  saturation  temperature,  °C  or  °F 

T*aii  Test  section  outside  wall  temperature,  °C  or  °F 

Urss  Root  sum  of  squares  uncertainty 

V  Voltage,  volts 

Vjjop  Test  section  voltage  drop,  volts 

pe  Electrical  resistivity,  ohm-m  or  ohm-ft 
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